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Abstract. "Rubber" coated rolling bodies satisfy a no-twist in addition to the no slip satisfied by "marble" coated bodies 1281 . 
Rubber rolling has an interesting differential geometric appeal because the geodesic curvatures of the curves on the surfaces at 
the corresponding points are equal. The associated distribution in the 5 dimensional configuration space has 2-3-5 growth (these 
distributions were first studied by Cartan; he showed that the maximal symmetries occurs for rubber rolling of spheres with 3: 1 
diameters ratio and materialize the exceptional group Gi). The 2-3-5 nonholonomic geometries are classified in a companion 
paper |29'| via Cartan's equivalence method [ 19 1. Rubber rolling of a convex body over a sphere defines a generalized Chaplygin 
system |21 57 , 35 27 39 1 with 50(3) symmetry group, total space Q = 50(3) x 5 2 and base 5 2 , that can be reduced to an 
almost Hamiltonian system in T*S 2 with a non-closed 2-form a^H- In this paper we present some basic results on the sphere- 
sphere problem: a dynamically asymmetric but balanced sphere of radius b (unequal moments of inertia /,• but with center of 
gravity at the geometric center), rubber rolling over another sphere of radius a. In this example g>nh is conformally symplectic 
1591 : the reduced system becomes Hamiltonian after a coordinate dependent change of time. In particular there is an invariant 
measure, whose density is the determinant of the reduced Legendre transform, to the power p = — 1). Using sphero-conical 
coordinates we verify the results by Borisov and Mamaev 1 1 211 1 31 that the system is integrable for p = — 1 /2 (ball over a plane) 
and p = —3/2 (rolling ball with twice the radius of a fixed internal ball). 



*On a CAPES-Fulbright visit to Caltech, Winter 2005-2006. 
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1 Introduction 

Given a riemannian manifold Q" and as<n dimensional distribution D of subspaces of TQ, two different theories apply. One 
is subriemannian geometry [51 1 which appears in the study of underactuated control systems and gauge theories. Nonholonomic 
geometry (2l l24ll53l [5l l25l follows d' Alembert's principle, and describes mechanical systems with nonintegrable constraint^. 

This an a companion 1291 paper belong to nonholonomic geometry. We use the following terminology: "Rubber" coated 
bodies means adding the no-twist condition to the usual no-slip constraints of "marble" coated bodies. The dynamics for the 
former lives on a 7 dimensional phase space, 8 dimensional for the latter. "Waxed" marble bodies have no velocity constraints: 
skiding dynamics lives on a 10 dimensional phase space. 

Strangelly enough, while marble bodies have been extensively studied in the NH-literature, rubber rolling seems to have 
been neglected. This is curious because they should be easier to study and furthermore, rubber rolling has an appealing differ- 
ential geometric interpretation. 

In this paper we provide details of some results on rubber rolling over a sphere announced in [28]. Details about the 
classification of 2-3-5 nonholonomic geometries via Cartan's equivalence [19j will be submitted elsewhere. The expert can 
go directly to section 3, or pass quickly over section 2. Although we tried to be reasonably self-contained, we use freely the 
geometric mechanics jargon, see eg. |T1|2"1 |49||54"1 . We also presented an extensive but focused list of reference^ 

Using sphero-conical coordinates we verify the results of Borisov and Mamaev 1121 1131 that the system is integrable for 
b/a = and b/a = —2. A table of results similar to that of Borisov and Mamaev (tables 1 and 2 of 1 10|) for marble bodies 
rolling over a plane or a sphere) is in order. 

1.1 Extrinsic differential geometry of rubber rolling 

In this section we discuss the implications of rolling a surface £2 without sliping or twisting over a surface £1 from the extrinsic 
differential geometric perspective. In [29 | we show that this process belongs actually to their intrinsic geometries. 

Given a curve C in a surface Z, define the (C,S) adapted frame as the moving frame F = (t,u,N), where t is the tangent 
vector, A' the normal to the surface (so an orientation is chosen) and u = N X t is the surface normal to the curve. Denoting 
' = d/ds the derivative with respect to arc length, as it is well known we have 

t' = K g u + K„N , u' = —Kgt — x g N, N' = — K„f + T g u (1.1) 

where K g = t' ■ 11 is called the geodesic curvature, K„ = t 1 -N is the normal curvature, and X g = N r ■ u the geodesic torsion. 
Recall an elementary result from classical differential geometry 

Proposition 1 (see e.g. Struik, i58V . p.201, exercise 2, section 4-8.) The geodesic curvature of a curve C on a surface S is equal 
to the ordinary curvature of the plane curve into which C is deformed when the developable surface enveloped by the tangent 
planes to S along C is rolled out on a plane. 

Remark 1 The inverse problem is hard: finding a curve with prescribed geodesic curvature on a surface S with metric ds 2 = 
Edu 2 + 2Fdudv + Gdv 2 gives rise to a nonlinear second order equation for u(s),v(s). 

Writing F as the orthonormal matrix with colums t,u,N, then F~ l F' = A, where A is the skew symmetric matrix 

/ -Kg -K„ \ 

A = K g ig 

\ K„ -X g / 

Let us write the structure equations dl.lt as 

t' = (£>f X t , U = COf X u, N' = COf X N , COf = -Xgt-K„U + K g N . (1.2) 

This yiels a nice dynamic interpretation: as the curve is traversed with unit velocity, the frame is rotating instantaneously around 
the vector Cfl/r with angular velocity \ \(£>f \\- 

Consider two surfaces £2 and Si an suppose that S2 moves always "touching" £1 . We neglect physical intersections that 
may occur if one of them is not convex). At the contact point the normals will be equal or opposite (depending on which 
orientations one choses). We use the following conventions: for a closed convex £2* the exterior normal defines its orientation; 
if Z] is also convex we take the opposite orientations when they roll externally to each other. 

1 Hertz 1 34 1 introduced the word nonholonomic to designate nonintegrable distributions and was the first to call attention that, although 
using the same ingredients, the two theories are quite different. 

2 We apologize for omissions, specially on nonholonomic reduction and special geometric structures (like almost Poisson and almost Dirac); 
we just mention some new connections of nonholonomic systems in robotics and control | 55 |,| 62 |. 
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When there is no sliding, the tangent vectors at the corresponding curves of points of contact have the same length, so we 
may use the arc length as a common parameter for both curves. Let C\ :x = x\ (s) and Ci : x = X2(s) denote the corresponding 
curves of points of contact in Zi and T^. an d let us describe the rolling action 

g(s) = (R(s),a(s))eSE(3) (1.3) 

(g acts on x £ R 3 by g ■ x = Rx + a, meaning first rotate, then translate"). 

Let Cflg the angular velocity in space of the rotation, and assume the condition co^ x xi + R~ l d = 0. Then R±2 = x\ and in 
fact the corresponding adapted frames are related by 

R(s)F 2 (s)=F 1 {s), (1.4) 

i.e, they match under the action of R(s). Here we orient £1 so that the normals to the surfaces point in the same direction. 
Differentiating, we get 

F 2 _ 1 {R- 1 R')F 2 +F2 1 F{ = Ff 1 F[ 

There is no loss in generality in assuming that at the point of contact corresponding to s = we have F\ (0) = Fi (0) = id (so 
R(0) — id, and x\ (0) = ^2(0) = 0. Therefore, if R'(0) is the skew-symmetric matrix B given by 

B\2 = -<*>3 , &13 = &>2 , "23 = -<»1 , 

then 

co Fl = co F2 + co fi 

which gives the following relation^ between the invariants of the curves C\ and C2 : 

Kg'' = Kg + CO3 , K^ 1 ' = 1Q, - C0 2 , Xg = Xg - COi (1.5) 

A very important consequence for our purposes is the following 

Proposition 2 Under the condition of no twisting (CO3 =0) the geodesic curvatures of the contact curves at the corresponding 
points are the same. 



1.2 Main results 

Rubber rolling of a convex body <B over a sphere of radius a yields a generalized Chaplygin system, see 1271 1311 l39l 1571 . 
We have a principal bundle with total space Q = 50(3) x S 2 and base space B = S 2 . The symmetry group G is 50(3), acting 
diagonally. The metric is left invariant, and the constraints define a connection on the principal bundle, called the "rubber 
connection". 

The theory developed in |27| implies that the dynamics reduces to T*S 2 , with a non-closed 2-form CO,,/, = &J C J' + (J,K) 
and a "compressed" Hamiltonian. The term (J, K) is a semi-basic form, where J is the momentum map and K is the curvature 
of the connection. The configuration variable in the base, ye 5 2 , has a nice geometric interpretation. It is the Poisson vector, 
namely, minus the external normal vector of the rolling body at the contact point (or equivalently, the normal vector to the base 
sphere or plane), seen in the body frame. 

If a solution y(t) is found, then the angular velocity £l(t) of ® in the body frame is the horizontal lift of y via the rubber 
connection. To reconstruct the attitude matrix one needs to solve a linear system with time dependent coefficients, R = R[£l(t)]. 
To complete the reconstruction, note that the contact point in 5 fl is q(t) = aRy(t). 

In principle, the curve q(t) in S a could also be reconstructed from Q(f) = dN~ l (— y(t)) on 5/, (where dN s is the Gauss 
map of 2) using the fact that the geodesies curvatures Kg of q(t) and Q{t) at corresponding points are the same. If one 
succeeds to solve this nonlinear equation, then the rotation matrix R can be computed, by noting that it takes the adapted frame 
(Q/\Q:Q x y/\Q x y!> ~ Y) along Q{f) to the corresponding adapted frame (q/\q,q/a x q/\q, q/a), (modulo a trivial change of 
signs depending on the chosen orientations). 

In passing, we observe that while it is simple to reconstruct a curve in the plane from its curvature k(s), for a general surface 
this task results on a complicated nonlinear equation (for K f = this is already the geodesic equations). For a spherical curve 

we get a linear system with variable coefficient^ 4r = ox ^r- — \q ■ 

In this paper we discuss the example where the surface of the rolling body « is a also a sphere, of radius b. It is a "Chaplygin 
sphere", meaning a sphere of mass p where the center of mass is the geometric center, but dynamically asymmetric, i.e., unequal 
moments of inertia Ij. The metric is given by 

2T=p(l±b/a) 2 \\q\\ 2 + (AQ.,n) , A = diag/ 1 ,/ 2 ,/ 3 (1.6) 



3 We learned this result from Mark Levi 1451 . 
4 Is this solvable by quadratures? 
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where q EE S a is the contact point and Q., CO are respectively the angular velocity of the rolling sphere Sb with respect to body and 
space frames, respectively. The constraints are given by 

(l±b/a)q = (±b/a)axq , q (£) = 0. (1.7) 

(the plus sign corresponds to the external case). 

Hence, when the rolling body is also a sphere, the constraints are invariant under the right action of 50(3) on the second 
factor of S 2 x 50(3). This system is (morally speaking) akin to a LR nonholonomic Chaplygin system 1311 . For a true LR system 
with 2 dimensional base space one can guarantee the existence of a function / : 5 2 — » R such that rf(/co„^) = 0. Moreover, the 
conformal factor / would be. f = F~ l l 2 , the inverse of the square root of the determinant of the reduced Legendre transform 
(identifying TS 2 = T*S 2 ), where 

F(y) = d e t(L^) = ^^ . (1.8) 

a y hh hh hh hhh ) 

Indeed, we will show that in the rubber sphere-sphere problem, the reduced system is Hamiltonizable. However, we get a 
different exponent: 

f a ,b(y) = F ^> y^ s2 - d- 9 ) 

Probably an explanation for this mysterious exponent will come from a study of Chaplygin systems of the form G °-> 
G x H /G — » H/G where the base space is a homogeneous space, and G acts diagonally in the total space. This theory should 
build up fom the LR systems H ■— > G — > H /G studied by Fedorov and Jovanovic [31]. 

Borisov and Mamaev [ 12] observed that by taking suitable combinations of the reduced variables, solutions of rubber rolling 
of a sphere over a plane can be mapped to the solutions of Veselova's system 16 II 1601 . They have also shown integrability 1131 
(in this volume) in the case b = —2a, where the fixed ball has half the radius and is internal to the rolling ball. We confirm 
their result by showing that the sphero-conical coordinates separates the Hamiltonian in the new time. Using sphero-conical 
coordinates/; < A,[ < I 2 < X% < h, given by 



(l 2 -h)(h-^2) (h-h)(I 3 -X 2 ) 



(1.10) 



{h-h)(h-h) ' {h-h){h-h) ' (h-h)(h-h) 

Quantitatively, our main result can be given as follows: 

Theorem 1 Let (q,R) £ S 2 x 50(3) the configuration space coordinates, where q is the contact point in the base sphere S a and 
R the attitude matrix of the moving sphere 5/,. Lety = R~ l (q/a) = —Q/b be the Poisson vector, where Q is the contact point in 
the moving sphere, seen in the body frame. The equations of motion for J are governed, in a new time scale X such that 

dx/dt = F{f, b -^ , F^) = (l + ^) 2 (Xi+pb 2 )(X 2 +pb 2 ) (1.11) 
by a Hamiltonian system (the nonholonomic vectorfield is X„f, = F(y)~ Xh) 

X H : ^=H P , *L =-H x , 2H= P l , + , (1.12) 

dX d% [(^+^2)^2+^2)]^ [(X 1+ pb 2 )(X 2 +pb 2 )]^C 2 



where 



Cl = h l+ if {X2 _ Xl) — — C2 = i (1+ ^ 2 _^) — i±±^l — ( i.i3) 

4 V a' y 11 Q*-h)(h-*.\)(h-h) 4 V a' V U (1 2 - h){\ 2 -/ 2 )(Z 3 -X 2 ) 



Corollary 1 The terms X 2 +pb 2 in c\ and X\ +pb 2 in c 2 are a nuisance for Hamilton-Jacobi separation, but they disappear in 
two cases, discovered first by Borisov and Mamaev [13], One is seen immediately, (b — a)/a = — 1, that is, b/a^ 0, the planar 
case. The other is (b — a) /a = —3, where we have the cross factors (\ 2 +mub 2 )/(X\ — X 2 ) in P 2 and (X\ +mub 2 )/(X\ — X 2 ) in 

r 2- 



P\; to see that that the Hamiltonian also separates in this case, multiplying both sides of tl,12i by 



X 2 —X\ 



(Xi +mub 2 )(X 2 -\-mub 2 ) X\+mub 2 X 2 +mub 2 

It is not known if there are hidden integrals for other values ofb/a and inertias Ij (except for the case of two equal inertias), 
but this is unlikely. The case b = a is special because the constraints are holonomic, so we have directly a two degrees of freedom 
Hamiltonian system. It was somewhat frustrating to realize that in this case the problem does not separate in sphero-conical 
coordinates and in fact seems to be chaotic from numerical experiment^ 



We admit losing a bet to Ivan Mamaev and Alexei Borisov, and we owe them a dinner in a barbecue house in Rio. 
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2 Preliminaries 

2.1 Skiding dynamics 

As we all know, the boundary of a strictly convex body <B is a closed surface £2 with strictly positive Gaussian curvature at all 
points. 

We say ® is in standard position when is the center of mass and the principal axis of inertia are aligned with OX,OY,OZ 
with moments of inertia I\ < I2 < I3, respectively. Under a Euclidian motion g = (x,R) e SE(3) the surface goes to g ■ £2, the 
center of mass of g ■ % is at x 6 R 3 and the three principal axis become, respectively, the columns e\ , e%, £3 of the attitude matrix 
R. 

The Gauss mappings Nj : £/ — * S 2 are important ingredients for the sequel. We make no special requirement on A'i , but 
assuming £2 to be convex is useful, it guarantees that N2 is one to one and onto 5 2 . 

The reader can easily sketch a figure to visualize our notation. Let £1 be fixed and £2 move under the action of SE(3). 
The unconstrained configuration space is ten dimensional, (SE(3) x £2) x Ej . When we impose the condition that the two 
bodies touch, the configuration space becomes five dimensional. One way to do the dimension count is to consider equations 
F\(xi,X2,Xi) = for £1 and ^(Gl , &, £b) = for £2- The touching manifold M 5 C (SE(3) x £2) x £1 consists of solutions 
for the 7 equations in 12 variables 

Fi(X) = 0, F 2 (Q) = 0,g.Q = X, and VF l {X)/\F l {X)\ = VF 2 {g- x Q) . 

We have 5 effective equations as the last three count as two. In the sequel we will identify M 5 = £1 x 50(3) , where the first 
component is the contact point, and the second component the attitude matrix of the rolling body. Given q\ 6 £j and R e 50(3) 
we can find the corresponding point O2 £ ^2 an d the position x = x^m °f the center of mass. Indeed, since the normals are 
aligned, we have A'i (q\) = —RNz(Qz) , where we orient £2 with the exterior normal. For visualization convenience, we orient 
£1 in such a way that at the contact points the normals point in opposite ways (for instance, in the sphere-sphere case, both are 
the exterior normals). Hence 

Q 2 =N 2 - 1 (-R- 1 N 1 (q 1 )) (2.1) 
x = x CM = 91 -RQl = 91 -R(N 2 - l (-R- l Ni (qi))) (2.2) 

(in the sequel we may drop the suffixes when no confusion may arise). 

As found long ago by Euler, the kinetic energy of a rigid body <3 is given by 

2T(x,R,x,R)=p\\x\\ 2 +Ii0. 2 1 +I 2 O 2 ,+I 3 0,l , (2.3) 

where p. is the total mass, Ij the moments of inertia about the principal axis ei,e2, e 3 attached at the center of mass x, and 
O. = (i2j ,0.2,0.3)' the angular velocity vector written the body frame. More precisely, R~ l R = [O] is the skew symmetric 
matrix 

/ -o 3 o 2 \ 
[n] = q 3 -fi] 

\ -0 2 Oi / 

See Arnold 1 1 1 for lower/upper case notations: The angular velocity in space is RR~ l = [00], with CO = RO. 

There is no coupling between the translational and the rotational motions if the body is moving freely in space. 

What is the dynamics when we impose the (holonomic) constraint of skiding? For simplicity, we ignore potential forces. In 
order to compute the Lagrangian 

T{q h R,q h R) on Q = £i x 50(3) (2.4) 

governing the holonomical system defined by the touching conditions 02.1 12. 2t . we replace in 02. 3t the velocity x of the center 
of mass by 

-* = 9l-^(R, 9l )[^(^ 1 (-^ 1 ^i(9i)))](^.'?l)- (2.5) 
This task is not so simple, because we need the derivative of the map 

F : {R, qi ) (E SO(3) ^ R(N 2 -\-R- l Ni( qi ))) (ER 3 

appearing in J2.5b . This notwithstanding, all the geometric information we need is contained in the Gauss maps: 

Lemma 1 

dF \(R.m) ■(^»9i)=<»x (91 -x) + (dN gTa )- i (a>xJVi(9i)-tflVi(9i)-9i) (2.6) 

Proof. This is a simple exercise on Advanced Calculus. Here the rotational velocity R is written in the space frame, RR~ = [co]. 
Denoting by g(t) -£2 the current (located) position of the moving body, g = (R.x) 6 5£(3), we can reinterpret some objects that 
appear in the derivation, for instance R~ l dN2(Q2)R = dN g z 2 (qi) . 
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Example: sphere-sphere skiding. Let Ej = 5 fl a sphere of radius a centered at the origin, £2 = Sb a sphere of radius rolling 
over Si . Then \2.2\ becomes as expected to x = ( 1 ± b/a) q\ The plus sign corresponds to the external case. The Lagrangian is 

IT =p{\±b/a) 1 \\q\\\ 1 +I l a]+I 2 ti L 2 +h£ll (2.7) 

and the dynamics uncouple. Solutions are great circles in the sphere and Euler rigid body motion for the moving ball. 

Proposition 3 For the problem of a convex body H with surface £2 skiding over a sphere £1 = S a there is 50(3) symmetry. We 
have a principal bundle structure with total space Q = S a x 50(3), 

50(3)^0^0/50(3) , R-{q,S) = {R(q),RS) (2.8) 

given by the diagonal action. 

Along the fiber (2 is been moved rigidly rotations around the center of sphere S a . In a sense motion along the fiber is "pure" 
skiding. The base of the bundle is the sphere S a and we can normalize to the unit sphere, parametrized by y. There is a 
global section, namely q\ = ay £ 5 2 1— » {ay, I) e S a x 50(3) . Notice that the induced global trivialization is not given just by 
the direct product structure. Fibers are diffeomorphic to 50(3) via R (aR(y),R). Along the global section, q\ = ay is the 
contact point in £1 = 5 fl , Q2 = (— Y) is tne corresponding contact point in £2 (in its standard position). Given a curve 
R(t)eSO(3),R{0)=I,R(0) = [p], 

K(P) = 4 R(t)(ay,I) = (apxy,p) (2.9) 
dt |f=o 

gives a vertical vector. 

Remark 2 Reduced equations for the skiding dynamics can be obtained by Marsden-Weinstein procedure H50V . For the 
symplectic reduction we need the momentum mapping J : T*Q —> R 3 . An element ofT*Q is of the form (pR,Pq,), where 
fHHm£ .vO(3)* = R 3 via right translation to the identity, and p qi £ T*(S a ) can be represented by a vector r perpendicular 
to q\. Using the abstract nonsense rule " J of p at q , on X = p at q on V(X) " we get 

J:T*(S a xSO(3))^sO(3)*(=R 3 ) , J(p R ,p q ) = m + qi xr = £ . (2.10) 

The six dimensional reduced symplectic manifolds are M 6 = J~ l (£)/$ , where 5 1 denotes the isotropy group of I. 

2.2 No slip and no-twist constraints: kinematical relations for the Poisson vector 

The Poisson vector 

y = R- l N l (q l ) = -N 2 (Q2) (2.11) 
is the normal vector to £] seen in the body frame and is the basic object for rubber rolling. It is related to Q2 by minus the Gauss 
map. We get immediatelly from RQ2 = q\ '■ 

Lemma 2 Kinematical relations for the Poisson vector, no slip: 

y+£lxy = R- 1 dNi(qi)-qi = [R- 1 dN 1 (q 1 )R]-Q2 ■ (2.12) 

In this formula, 171 = (iVi )~' (Ry) when the local inversion is possibl^E I n particular, if £1 is a sphere of radius a or a plane 
(a = 00), the term R~ 1 dNi {q\ )R is just 1 /a times the identity, the kinematical relation becomes 

y+£lxy = R- 1 (q l /a) = -Q 2 = —-d(N2)y i -y (2.13) 
a a ' 

Furthermore, if £2 is also a sphere, of radius b > 0, 

y+KQxy=0 , k=— (2.14) 
a±b 

where the plus sign corresponds to the external case, minus when one of the spheres contain the other. We can use only the 
plus sign and allow b to be negative to represent the internal case (if b < and \b\ > a then the fixed sphere is inside the rolling 
sphere). 



Points on E] where one of the principal curvatures vanish may be specially relevant for the dynamics. 
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2.3 Rubber connection and curvature 

Imposing the no-slip constraints produces a map 

<Di->?l (at qi ) (2.15) 
which gives a distribution of 3-subspaces in Q. To rule out twisting, one adds 

(a-N\(q\) =0 or equivalently Q. ■ y = . (2. 16) 

In more detail, the no-slip constraint in Q = £j x 50(3) follows by replacing the left-hand side x in ( 12.51 ) by 

x= (OX (x — qi) . (2.17) 

In view of lemmaQ] the result is 

A I J 



<h = dN'l (co x AT] ( 9J ) - dN\ (gi Hi) (2. 18) 



or, in a more symmetric form, 

[dN g z 2 +dNi](q 1 )-q 1 =(oxNi{qi) = mx(-RN 2 {Q2))- (2.19) 
Since for no-slip the tangent vectors of curves q\ (t) and Qi (t) correspond, RQi = 4\ , equation d2.19t gives 
Q x y = R- 1 [dNgZ, + m ] («l ) ■ 91 = dN 2 (Q 2 ) ■ & + ««i (91 ) • 4i ■ 
Imposing in addition to J2.19I ) the no-twist constraint ( 12.161 ) allows us to solve for co in terms of q\ ■ 
Lemma 3 The " rubber connection" : 

CO = Ni X ( [dN gL2 + dN x ] ( 9 i Hi . (2.20) 

This formula defines an Ehresmann connection on the bundle 50(3) X £1 — » £1 . Equation J2.20b can be rewritten in the body 
frame as 

n= [(Id+{R- l dNiR)(dN 2 r 1 )}(y)xy (2.21) 

Note that there is no 50(3) equivariance unless Ej is a sphere. 

Consider a small curve C\ around a point q\ £ £1 . The curvature of the Ehresmann connection is the limit 

lim ( how much a frame attached to £2 rotates ) / ( area inside the curve in £j ) , as C\ shrinks . 



We have computed the curvature for the case of a sphere rolling over another sphere, see section [331 below. In order to 
appreciate the difficulties of the calculation, we challenge the reader to attempt computing by brute force the curvature in the 
case of a surface of revolution rolling over a sphere. We claim, however, that it is possible to give a formula for the Ehresmann 
curvature of the rubber connection in the case of general surfaces £1 , £2 in terms of natural geometric objects. Accepting at face 
value the information that rubber rolling belongs to the intrinsic geometry 1 14], the idea is to use proposition|2]together with the 
Gauss-Bonnet theorem. The result (see details in the companion paper [29 1) is a kind of dynamical "Egregium theorem": 

Theorem 2 (Bryant and Hsu 114]). Curvature of Ehresmann connection, general case: Let kj{qi) the Gauss curvature o/£; at 
the corresponding points q,, such that q\ = gqx , g = {R,a) £ SE(3). Then 

K=(l-'^4)h(q l )<n: 1 . (2.22) 

where d"L\ is the area form of £] . 

This result is remarkable because the connection 1-form already involves the derivative of the Gauss maps, so at first sight the 
Ehresmann curvature should involve second derivatives. 



Remark 3 Koon and Marsden have shown that a key ingredient for the equations of motion for nonholonomic systems is the 
curvature of the local Ehresmann connection associated to a splitting of a coordinate system (see [6. 431). Both for marble or 
rubber rolling this splitting is global. See Theorem\3\below for the case when £1 is a sphere. 



2 PRELIMINARIES 



9 



2.4 Equations of motion: general physics approach 

In this section we follow Borisov and Mamaev 1101 closeljQ. In addition to the kinematical relation i2. 12t . three differential 
equations for the rolling body can be derived from 

(CM = *CM > 

where £cm > X CM are respectively the angular momentum and the torque of the external forces with respect to the center of mass. 

In basic mechanics textbooks it is sometimes overlooked that the above formula does not hold in general if one takes q\ 
as the base point instead of the center of mass CM. Recall that the angular momentum with respect to the contact point q\ and 
with respect to the center of mass x are related by 

t qi =(x-q 1 )x t a + l C M - (2-23) 

Differentiating, we get 

4 = [( x - 1\ ) x fext + ^cm] + -j- [(x- qi)} X fix , 
and we recognize the term (x — q\ ) x f ext + t^m as the torque about the contact point. The correct form of the torque equation is 

4 =T 9 , + —[{x-qi)]Xfd. (2.24) 

containing an ammended term, which can be called the dynamical torque, 

d 

Idyn = jdC* _ 9l)J X/« . (2.25) 
The non-sliping condition implies x = CO x (x — q\ ), so in the space frame 

4 = % + — [(x- qi )] X fia X (x - q{) . (2.26) 

If we go to the body frame, we ge0 
where we recall for clarity, that 

Eqi = RL, x - qi = RQi, co = RCl, z = RT . 

In order to obtain a linear relation between L and £1 (depending on the Poisson vector y), we rewrite the total energy of the 
system as 

2H = 2T = {l qx , co) = (L, Q.) , whith RL = l qx , 

where t qi is the angular momentum with respect to the contact point viewed in space, L its coordinates in the body frame. 
Substituting the no-slipping constraint in l !2.3t , one gets 

L=A£l + uQx(axQ) , A = diag{I u I 2 ,h) ■ (2-28) 

L = , is the angular momentum of the rolling body with respect to the contact point, seen in the body frame. Relation J2.28t 
can be inverted, 

£1= (A+v\\Q\\ 2 id)- 1 L + a^QjiA+vWQlfid)- 1 (Q) (2.29) 

where 

(6,(A+A/iieii 2 ^r^) 

'\-n{Q,{A+n\\Q\\ 2 id)- l Q) 
It will be convenient to use the shorthand A = A \Q\ \ 2 id. 



L + axL+(axQ 2 )x{ijQ2) = T (2.27) 



Proposition 4 (Counting dimensions and equations) In the no-slip case we have an eight dimensional phase-space. Using 
(R,y,L) £ SO(3) x S 2 x R 3 as coordinates (we are assuming that q\ = (N\)~ l (Ry) can be locally inverted), the dynamics is 
given by the momentum equation \2.27\ , the kinematic relations \2.12\ and the attitude equation R = R[Q]. We will have a 
closed system of equations since fl is a function of (L,y) via formula \2.29\ . If we add the no-twist constraint, the phase space 
is seven dimensional: one adds the constraint (£2,y) = 0, and uses it to eliminate the unknown torque around the normal. 

Reconstruction of q\ motion in Zi can be done either via q\ = RQi, (where as we recall, Q is related to y via the Gauss 
map) or simply inverting N\ (q\ ) = Ry. In general, the kinematical relation couples R with the other variables, but when is a 
plane or a sphere, the attitude equation decouples from the equation for y and L, a consequence of 50(3) symmetry. 



7 Formore abstract approaches, see eg., |42|, |46],|47|. 
8 This is precisely eqs. (1.1) of 1101 in our notation. 
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Remark 4 Ifx and q\ move in such a way that a combination of them with fixed coefficients remains constant, then x and q\ are 
proportional. Consequently, the ammended term in \2.25\ vanishes identically. For instance, this happens in the case of spheres 
rolling over spheres (or a planef^ 

Sphere-sphere rubber rolling. As the dynamical torque is identically zercP^I , we are left with 

L + £lxL = zy, y=Kyx£1 , L = ASl+pb 2 (Sl,y)y = AQ. , with (n,y)=0 (2.31) 

with K = a/(a + b). The right hand side T is the external torque about the contact point in the body frame (zero for the no-slip 
case, but twisting allowed.) In view of the no-twist constraint (£2, y) = it is immediate that Q. = A~ 1 L, and we get the following 
expression for the multiplier: 

t(Y,A _i y) = (AT 1 ! xL,l~ l i) + K{AT l L,l~ l Lx y) . (2.32) 

The problem of a rubber ball rolling on a plane corresponds to K = 1 and was discussed briefly in [27 ] . We showed the 
existence of an invariant measure 

v = F(yy l/2 dL l dL 2 dLj,djidy2dj 3 (2.33) 

with 

F=W3f(Ar)^^ + trR^) • (2.34) 

\ hh hh hh hhh J 

Borisov and Mamaev 1 12 1 called the attention that the solutions of the rubber ball on the plane can be set into correspondence 
with solutions of Veselova's system |61 , 60 1 namely, a rigid body with one right invariant constraint (formally, b = 0). They 
used a trick that seems to go back to Chaplygin, which is really a "Columbus egg": one takes suitable linear combinations of y 
and L and show that the equations of motion correspond. Unfortunately, this clever maneuver only seems to work when K = 1 . 

2.5 Marble Chaplygin sphere over a sphere or a plane 

This section is intended as background material, just for comparison with the rubber rolling problem discussed next. For more 
details on what is known about marble rolling of convex bodies over a plane or a sphere, see 1 10 1. 

In 1903 Chaplygin [22] integrated using hyperelliptic functions the following system for (L,y) 6 x 9? 3 : 

L + £2xL = 0, (2.35) 

y+flxy = 0, (2.36) 

where SI is defined by jl.29\ with Q = rj, namely, 

n(L,y) = (A+ pr 2 id)- 1 L + a(L,i)(A+ pr 2 id)- x (y) . (2.37) 

with 



1 — iir 2 (y, (A + ixr 2 id)" 1 y) 

They describe the motion of a marble ball of radius r and mass p., moments of inertia Ij about the geometric center (which is 
also the center of mass), rolling without sliping over a plane. Twisting motions (rotations about the vertical) are allowed. L is 
the angular momentum of the ball with respect to the contact point, SI the angular velocity and y the vertical vector, all viewed 
on a reference frame attached to the ball. 

This system has four independent integrals 

h = (L, Ot) , h = (L,L),ft = (y, y) , f 4 = (L, y) (2.39) 

where 2/i is the energy and f\, called the area integral, is the third spacial component of the angular momentum. Topologically, 

the common level sets of these integrals are tori; Chaplygin showed that there is a smooth invariant measure 

with 

F(j) = (Y,A _1 y) > A=A+pr 2 id (2.40) 

pr L 

Integrability follows from the celebrated Jacobi 's last multiplier method, which says that an ODE on a two dimensional manifold 
having a smooth invariant measure can be solved by quadratures. Detailed remakes of Chaplygin's paper appeared recently, see 
Kilin 1 37 1 and Duistermaat 1 26]. The planar motion can be reconstructed by integrating x = r(j)(t) x k. 



'Historical note: Chaplygin seems to have been the first one to discuss geometrical conditions under which the second term vanishes. See 
I53| |21 I4T1 for discussions on this often neglected issue. 

l0 For abstract explanations, see discussions on nonholonomic Noether theorem in (2 
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Remark 5 The issue whether Chaplygin marble system is hamiltonizable or not has been discussed in fiE\ [9] \26\l . Perhaps 
due to some misprints we could not verify the Borisov-Mamaev bracket, but we are glad to know that it has been verified and 
explained independently by Naranio R52]l . Theoretical issues about almost Poisson brackets are outside the scope of our paper, 
but we believe will be much in evidence in the next years, see [15. 16 47. 48 1 for background. 

To extend Chaplygin's sphere-plane equations to the sphere-sphere case, it suficies to replace equation d2,36t by the kine- 
matical relation 

Y+Kflxy = (2.41) 

with K = a /(a + b) (Recall that b > in the external case, so < K < 1 ; in the internal case where b < and \b\ < a, K > 1 ; 
and K = 1 for rolling over a plane (a — » °°). One can also consider the case — °° < k < 0, when the fixed ball S a lives inside the 
rolling ball S\ b \,b<0,\b\ > a.) 

The formula for the energy is the same as in the sphere-plane case, 

2// = 2r=^ 2 ||flxY|| 2 +/,flf +/ 2 fl2+ / 3 fi 3 • (2-42) 

We now give an argument to explain why the invariant measure stays the same, independently of the value of K. 

Proposition 5 There is a smooth invariant measure for the marble rolling of a sphere over a sphere. The density function is the 
same as in the case of Chaplygin's sphere rolling over a plane, F(y) — with F given by i2.40i . 

Proof. As a warm up, let us verify first that the total energy function ( |2.42t of the sphere rolling over a sphere (a formula in 
which the radius of the base sphere does not appear) is conserved along the trajectories of J2.35ll2"!4~0 . This direct check will be 
instructive, since all sphere-sphere systems have the same energy function and the same L equation, but a family of y equations, 
parametrized by K. We still get H = grad^ „\ H ■ (L, y) =0 because both 

V L HL = 0, (2.43) 

V y H-y = 0. (2.44) 

Equation J2.43t is immediate, since H = fl and L = Q, x L is perpendicular to £1. Now, we know that the proposition holds 
for the case K = 1 . Equation J2.44t must be therefore true for K = 1 . Now, for an arbitrary K, the left hand side of l |2.44t only 
picks this multiplicative factor so it must also be zero and there is no need to do the computation). 

In order to verify that the measure is invariant, we follow a similar procedure. Let us scrutinize the derivation presented 
Duistermaat, | 26| (lemma 7.1) for the case K= 1. His approach was to split 

div (7X) (yx n(y,L) ,l x a{y,L)) 

as the trace of the derivative with respect to y of the first component plus the trace of the derivative with respect to L of the 
second component. The latter trace is obviously identically zero, whereas the first one is found to be (see (7.2) in that paper) 

div y ( Y x Q(y,L))) = ii j± '- . (2.45) 

On the other hand, the directional derivative of F(y) is found to be 

(gradF,Yxfl(Y,L)) = 2(yxA _1 y^ -1 ^) (2-46) 

and hence 

div (yX) [F-'/ 2 (yx Q.(y,L) ,L x Q(y,L))] = divy[F~ l / 2 yx £l(y,L) = . 
Now, comes the argument: the only change caused by the presence ofK is to multiply both (12.451 and \2.46\ by this same fator, 
so the net result stays zero. 

Functions f\,f2,fo are still integrals of motion for any K, but the "area integral" is lost! It is easy to see that it holds only 
for K = 1. In the three dimensional manifolds M 3 defined by level sets of the three (surviving) integrals we have an invariant 
measure. Flows on three dimensional compact manifolds with an invariant measure are of interest for dynamicists and one 
expect that the motion will be ergodic in M 3 except perhaps for some special values of k, where some new integral of motion 
may appear. 

Borisov found the only known (perhaps unique) new integrable case, K = — 1, corresponding to an internal fixed ball of half 
the radius of the rolling ball, see 1 10 1, pag. 194. 

Remark 6 The external case with equal radii (corresponding to K = 1/2) could be (erroneously) thought to be also integrable, 
for a simple reason. As we will see below, the two dimensional rubber distribution is holonomic for two spheres of equal radius. 
One could hope that relaxing the no-twist constraint would keep the distribution holonomic, as it amounts to allow any finite 
rotation of the moving sphere about the normal at any given (but arbitrary) point on the base sphere. Each leaf f would then 
be three dimensional, forming a a S 1 bundle over S 2 , and the topology would be (probably) that of SU (2). Unfortunately there 
is a serious blunder is this reasoning. Although span(Xi ,X2) is integrable, adding the infinitesimal rotations about the contact 
points produces a 3 dimensional distribution span {X\ , X2 , X3 ) with growth 3-5. 
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3 Rubber rolling over a sphere as a 50(3) Chaplygin system. 

3.1 Reduction to T*S 2 . 

Consider the principal bundle with total space Q = S a x 50(3), 

50(3)^2^2/50(3) , R-{q,S) = {R(q),RS) 

given by the diagonal action. It is not hard to see that the no-slip and no-twist constraints define a principal bundle connection 
for any rolling surface £2. We call it the rubber connection. 

Recall J2.20K the Ehresmann connection for the rubber rolling of a surface Z2 on an arbitrary base £1 : 

00 = Nt X ( [dN gLl + dN 1 ](q 1 )-q 1 . 

When is a sphere there is an extra feature, namely, equivariance, so we have a principal bundle connection. We need to 
show that if CO solves both equations i2. 16ll2~T9l l. which in this case become 

CO X qi = [/ + adN g z 2 (q\)]q\ , CO • q\ = , 

then R(£> solves 

RwxRq! = [I + adN Rg z 2 (Rqi)](Rqi) , R(O-Rqi=0. 
The second fact is obvious, and the first follows frorrF*! 

dN RgL2 (R qi ) = RdN^iq^R- 1 . 

See 1 39 1 for reduction written in Lagrangian form. The reduction procedure in Hamiltonian form, following the recipe given 
in 127] is as follows: 

Theorem 3 Let £2 = £ be arbitrary but £[ a sphere of radius a. Then rubber rolling is described by a non-abelian Chaplygin 
system on the bundle Q = 50(3) x 5 2 (fl), with symmetry group 50(3) acting diagonally. The constraints are given by the rubber 
connection 

(0=^-x[dN gT .- qi + ^-]. (3.1) 
a a 

Define the projection is by n(R,q) — RT^qja 6 5 2 , R € 50(3) , q £ S 2 (a). Then the dynamics reduces to an almost-symplectic 
system (H,G)nh) in T*S 2 . The Hamiltonian H in T*S 2 is the Legendre transform of the compressed Lagrangian T(y,y) given 
by 

2T = {AQ.,Q) with a = [{Id+-(dN 2 r l (y))}(j)xy , A = A+fi\\Q\\ 2 id , Q = -N L l (j) (3.2) 

a 

obtained from the original one by horizontal lifting via the connection. The 2-form is given by co„/, = C0^ fl) ^ + (J,K) and in 
general is non closed. The ammended term (J,K) is semibasic; J is the momentum \2.10\ of the 50(3) action and K is the 
curvature \2.22\ of the rubber connection. 

The domains of J and K are matched by the Legendre transformation from TQ to T*Q. Their Ad and Ad* ambiguities 
cancel out automatically. In some more detail, 

(J,K) = (l-^-)(y,M)dMeas>, (3.3) 

where y = q/a is the contact point in the base sphere, kj_(y) is the Gauss curvature of the rolling surface at the current contact 
point 2 = — Ny 1 (y), and M = AO. is the angular momentum at the center of mass of the rolling body. 

We plan to use this result as a basis for future work. The first in line is the rubber rolling over a sphere of surfaces of 
revolution with two equal inertias. For the rest of this paper we restrict ourselves to the sphere-sphere example. 

3.2 The rubber sphere-sphere connection 1-form. 

We show that for the sphere-sphere problem there is an invariant measure, and even better, the reduced system is conformally 
symplecticT 2 !. The expert can skip many of the redundant calculations below, intended as checks of the general results presented 
earlier. When S2 ls a sphere of radius b, d3.lt becomes 

co=^(i±~)<nx?i. (3.4) 

a L b 

Recall that for simplicity we drop the ± sign, using the convention that b is negative one of the spheres contains the other. 

"Latexing is more difficult than the result. A drawing should make it obvious. 
12 We conjecture that this is also true for bodies of revolution. 
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Lemma 4 The horizontal lift of y e T y S 2 at (aR(y),R) <EQ is given by 



j^Hor(j) = (-bR(j), (1 + -) (L«)*[tX Y]) ■ (3-5) 



where (L«)* indicates the left translation and [,] f/ie standard isomorphism from R 3 to sO(3). 
Proof. The formula looks odd due to the fact that the 50(3) action is diagonal. Note first that 

(n t ) qM) (q,R) = ^-(R-'q-R-'RR-'q) . 

TakeoK(y) = q,q = -bR(y) and R = (L R )*[y x y]. Then 

= -^Y-(1 + ^)(?xy)xy)=Y 

and a short calculation shows that l !3.4t holds: 

-L(l + ^)R- l (q ]X qi) = -(l + j)yx(-bj) = (l + -)yxy. 
a b a b a 

In order to find the connection form, it sufficies to compute along the global section yn (ay, I). 

Proposition 6 The connection 1- form (]) interpreted as a 1-form in 5 2 x 50(3) with values in sO(3) = R 3 , is given at the point 
(y, /) of the global section by 

<b(y,o) = (l+b/a)yxy-(b/a)[o-(l+a/b)(o-y)y] . (3.6) 
Proof. By a direct check one verifies, as desired, that 

<KPXYP) = P , <t>(--)Y,(l + ^)YXY) =0. 

While the metric is left invariant under the diagonal action of G = 50(3) on S a x 50(3), looking at J3.4t we observe that 
in the sphere-sphere problem the constraints are invariant under another action of 50(3) on5 2 x 50(3), namely the right action 
on the group component: (y,R)S = (y,RS). This makes the problem reminiscent to a LR nonholonomic Chaplygin system 
H ■— > G — > H/G, see 1311 , except that here the bundle structure is of the form 

G^GxH/G^H/G (3.7) 

with H = 5 1 the rotations about the vertical axis, H/G = S 2 the homogeneous coset space. 

Assuming that results similar to those for LR systems would still hold, since the reduced dynamics has two degrees of 
freedom, we predict the existence of a function / : 5 2 — > R such that d(fd}„i,) = 0. In other words, the reduced system in T*S 2 
would be Hamiltonizable, since the nonholonomic 2-form is conformally svmplecticP 3 ! and in particular, there will be a smooth 
invariant measure. Actually, we guess / from the density function of the invariant measure, relative to the Liouville volume 
of T*5 2 . The good news is that it is possible to verify directly if a nonholonomic system has a smooth invariant measure, see 

SEIEHIISIISI- 

In future work we plan to pursue a LR theory for principal bundles of the type d3.7l l. Here, we use an alternative, albeit 
tentative approach: in order to make the LR analogy more plausible, we consider the whole "onion" of base spheres S a ,a > 
(the singular limit a = represents Veselova system). Then our mechanical system lives on an extended configuration space, 
Q = E(3)' = 50(3) x (R 3 — 0) with a metric which is manifestly invariant under the left action of 50(3) : 

2T=p.{l±b/\q\) 2 \\q\\ 2 +l l Q. 2 +I 2 Q.l+I 3 Q.i 

Note that the presence of \q\ in the denominator forces us to remove the origin of R 3 . 
The pseudo LR system has bundle structurj^l 

50(3)) E{3)) -+ R 3 , R-(S,q) = (RS,Rq) (3.8) 

We now rewrite the no-slip, no-twist conditions in terms of differential forms. The right and left invariant forms in 50(3) 
will be denoted here 

/ -p 3 p 2 \ 

dRR' 1 = -pi ,R^dR^\ -A,] I . (3.9) 





13 See 1 32 33 59 63 1 to appreciate the consequences of being conformally symplectic. 

14 Using a more abstract language helps to generalize to configuration spaces that are semi-direct products of Lie groups. 
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The distribution £> of 2-dimensional subspaces in £(3)* with \q\ ^ 0, is the joint kernel of the four differential forms 
r|i, 1 < i < 4 in £(3) defined by 

(Tll,Tl2» T l3) f = P-rvi( 1 +u)<l xd <l J ^4 = q\dq\+q 2 dq 2 + ?3<*?3 • (3-10) 

where p = (pi,p2,P3)'. Notice that the Pffafian equation (q.p) = defining the no-twist condition is satisfied automatically 
when r|j = 0, 1 < i < 4. Notice also that 5 G 50(3) acting on the right on an element (R,q) £ SE(3) gives (RS,q). This is best 
seen using the costumary representation for E(3): 

<*.'>«( o 

This implies 

Proposition 7 Consider the distribution 2> ;'n £(3)* defined by the constraints T],- = 0, 1 < ; < 4. These forms are invariant 
under the right action of SO (3) in SE(3)', that is D(R.q) = ^(I.q) m R- 



3.3 Curvature of the sphere-sphere connection 

This section could be avoided in view of ( 12.22b but we decided to include it for completeness. We already know the connection 
1-form (]) at the global section. In order to compute the curvature 2-form we need an expression everywhere. This can be done 
by equivariance, but this is a somewhat unpleasant task [ §(q t R)(q,R) is obtained from i3.6i by computing R of the result of 
(j)(y, <7) with the inputs y «— R~ l q/a, y <— R~ l q/a — CO x y, [a] <— R~ l R, CO = RR~ l . Also, we must extend it to a connection in 
R* x 50(3) <-» SE(3) m — » 5 2 , by adding a component (q,q/\q\) to ()), and replacing a by \q\ and q by q — (l/\q\ 2 )(q,q). ] 

After putting everything inside (]), we would get an unpleasant formula. Computing the curvature directly would be a "tour 
de force". We prefer to guess the curvature by a wrong but algebraically elenant calculation. We then reinterpret and confirm 
the guess using elementary geometry. 

We can use cartesian coordinates in R 3 , since we extended the base space to an "onion" of spheres. Consider the vectorfields 
in 5£(3)* given by 

Y x = d/dx+ (- + I)(£* 2 - y -X 3 ) , Y 2 = d/dy+(- + h(--Xi + -X 3 ) 
a b a a a b a a 

where X\,X2,X 3 are the right invariant vectorfields of 50(3) ([Xi,X2] = — X 3 , etc). These vectors are horizontal at the point 

c la = (0,0, a), R = I. A simple mental calculation gives us 

2-(- + 7 ;)-(- + t) 2 ] = --t 2 - • 
a a b a b a^ o z 

Now, both Y\ and Yj_ must be multiplied by — b in order to be the horizontal lifts at q a of the unit vectors (1,0,0) and (0,0, 1) at 

the north pole of 5 2 = 5 2 (1), see equation d3 - 5 1 > in lemma[4] As a result, we have to multiply that result by b 2 . From H.Cartan's 

magic formula, we get (note the extra - sign): 

Proposition 8 The curvature of the connection, interpreted as a 2 form in the unit sphere 5 2 with values in R 3 = ,sO(3), is 

b 2 

*(Y) = (l-^)Yrfarea s2 (3.11) 
a L 

This calculation can be made more rigorous. We forgot to compute Y\§{Y-i) — l^'K^l) m H. Cartan's "magic formula", at 
(0,0, a). By a symmetry argument, one should be able to show that it must vanish. 

We now present a geometric way to compute the curvature. Consider a small circle C\ in the base sphere S a . Roll along 
C\, without slip nor twist, a ball 5;, (for concreteness, we do it externally), until it comes back to the initial point in C\. The 
curvature is the following limit: 

lim ( how much a frame attached to Si, has rotated ) / ( area of the cap in S a ) , as S a shrinks . 

Now, Proposition [2] tells that the corresponding curves have the same geodesic curvatures, so the contact point on 5^ must 
move in another circle After one turn arond C\, the end point in C 2 is not the same as the initial point. For instance, the 
marker does several turns plus a fraction if b is smaller than a. Moreover, proposition Q] allow us to replace the spheres by their 
enveloping cone 



15 Cone rolling has many applications, ranging from transmission gears to leisure stuff, see http://www.kurt.com/gear.html , (no 
relation with one of the authors), http://www.lpercent.com/store/cart/RLZENCONE.html (no relation with any of the authors, and 
http://kmoddl.library.cornell.edu/model.php?m=reuleaux (yes, we enjoy anything related to Cartan's 1:3 rolling! 
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The equivalent problem is the following: let a cone Z2 roll around another cone Zj until its generator returns to the same 
place. How much the position of a marker point Qi in S2 nas rotated around axis li ? We call this angle divided by 271 
the holonomy. The reader should be able to sketch a figure. The moving cone S2 r °H s around the fixed cone £j, having a 
common a different generatrix g at every instant. Draw right triangles VO\q and VOiq where V is the common vertex, 0, in 
the axis of Z; are the centers of the circles Q, and q is the contact point. Let a = LO\ Vq and (3 = LOtV q the cones apertures, 
6 = 7l/2 - a = LOiqV, = Jt/2 - (3 = LO q qV . 

The answer is easy: 

I holonomy I = length of the base circle C\ / length of Ci = sin(oc)/sin(P) = cos(0)/cos(<|>) . 
Now by elementary trigonometry, 

atan(e) = fctan(0) (= £) . (3.12) 
which by the way says that the geodesic curvatures are 

Kc,=Kc 2 = l/£. (3.13) 

Thus 

a sin(6) 

holonomy = -__L2 , (3.14) 
b sin((|)) 



Moreover, from ( 13 . 1 2b we can solve for 



sin((|>) = - 



tan(9) 



'l + f;tan 2 (9) 

so that 



/ a 2 

holonomy = -cos(0)\/ l + ^j tan 2 (6) . (3.15) 

Why the minus sign? The angle between the normals is a + P, and since + a = 7t/2, + (3 = 7c/2, both a and P tend to 7t/2 
when 6 (and so 0) tends to zero. So the two normals tend to become opposite to each other ! Actually, here is one more twist: 
a frame attached to the sphere S/, rotates by 271 — 271 1 holonomy |. The extra 271 is due to the "Phileas Fogg effect" (known to 
mathematicians as Hopf's umlaufsatz). Now, the area of the cap in sphere S a is 27Ta 2 (l — cos(6)) ~ 7ta 2 6 2 . 
In short, we end up needing to compute the limit 



2 1- Jl + fitan 2 (e)cos(e) 

1™ — T~ XT 

6^0 a 2 6 2 

This is a nice calculus exercise, the result being that the curvature is inded -\ — -X times the infinitesimal rotation about the 
normal around the contact point q in sphere S a . 

Remark 7 Among the ratios a/b, two are exceptional: a — b provides a holonomic distribution since the curvature vanishes. 
Carton has used the 3:1 ratio to construct his exceptional Lie group G2, see M8\ \7v. We hoped that the corresponding rolling 
problems would have some special features, but that was with no avail. 

3.4 Hamiltonization of the sphere-sphere system 

We will compute step by step the 2-form con// = dpydy+ (J,K), and the reduced energy function H : T*S 2 — > R, to convey the 
algorithmical nature of the method. 

• The momentum mapping is given by J(m q ,r) = m + q x r , see equation l |2.10| l. 

• Pairing the momentum with the curvature < |3 . 1 U along the global section q £ S 2 (a) 1— » (q,I), q = ay we get 

b 2 

(J,K) = (1 - - T )(Y.M)rfarea S 2 , M=A£l 
a z 

In order to match the domains, we follow the "clockwise diagram" presented in 1271 . This diagram is a consequence of a 
Hamiltonian counterpart of Lagrangian reduction |20 | using quasi-coordinates 1401 . 

• The Legendre transform TQ 1— » T* Q induced by the Lagrangian r = ^[^(l + ^) 2 |4l | 2 + {AQ.,Q.)], at the point (q\,R)^ i «) 
(v,£2)( flY ,/) e T (ayJ)Q' ls ( p v, p i)( qi ,R)> given by 

P I = M = Aa, P v = r = (i(l+ -) 2 v . 

a 
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• The Legendre transform of the horizontal lift Hor(j) = (v,fl) = (-by, (1 + |)yx y): 

m= (l + -)A(f xy) , r = ^(-fo)(l + -) 2 y 
a a 

The end result is 

(7,^) = (l-^)(l + -)(y,A(yxy)d5 2 = (l--)(l + -) 2 (y,A(yxy) fl !5 2 . (3.16) 
a z a a a 

• To obtain the compressed Lagrangian T re d, we compute the Lagrangian at Hor(j): 

^/ = i(l + -) 2 [^¥ + (A(yxy),yxy)] (3.17) 

• The Legendre transform in the base, py = dT rei j /3y = G(y) ■ y: 

Py = G(y) ■ y = ( 1 + b - ) 2 [A* 2 y + A(y x y) x y] (3.18) 

• The determinant of G(y), restricted to TyS 2 is 

F(j) = [dctLe gcomp ] = ( h I 2 I 3 )(l + b -) 2 (Vy,Y) +^[^ + ^ + + . (3.19) 

Only this last assertion represents a nontrivial calculation, but that was done in |27 1, for the limit a — > °° of the rubber ball rolling 
over a plane. In that paper, we announced ((3.38), Proposition 3.6 of 1271 ): 

Proposition 9 For a = °° ( rubber ball rolling over a plane ) 

d{f(a T cat +(/,*))= , / = F-'/ 2 . (3.20) 

Proof. This was checked directly by computer algebra, using spherical coordinates. It was useful to pull back the canonical two 
form in T*S 2 to TS 2 via py = G(y) • y. We can provide a Mathematica notebook for those interested. 

We now show how to adapt this result for b/a ^ 0. First note that the conformality condition can we rewriten, when the 
base space of the bundle has dimension two, as 

d{logf)A(O ca „+d(J,K)=0. 

From J3.16t and J3.18I ) we observe that (i) can , when pulled back to TS 2 , picks the factor (1 +b/a) 2 , whereas (J,K) picks 
(1 — b/a) (I +b/a) 2 , that is, one fator 1 — b/a surplus. Thus we can't have the same / as in the case a = °°, but if we change / 
to f x ~ b l a , that should work. 

Theorem 4 The reduced rubber ball rolling over a ball is conformally symplectic, 

d\f a ,b(y)(®caf + (J,K))}=0 (3-21) 

with conformed factor 

fa,b{i)=F b S ■ 

Proof. For the suspicious reader, we can provide a Mathematica notebook using spherical coordinates. But a derivation (done 
by hand) follows from sphero-conical coordinates, as we show in the next section. 

Remark 8 The presence of the exponent is somewhat mysterious, but it makes sense: when b = a, then f = 1 (as it should 
be) since in that case (J,K) = and the problem becomes Hamiltonian. The predicted exponent from the results by Veselov 
and Veselova H61V or Fedorov and Jovanovic [31 J would be — 1/2, but this corresponds only to b/a — » 0. We believe that the 
problem resides on the singularity at the origin a = 0, and the mystery will be resolved with a thorough study of bundles H3. 71 ). 

All these problems have symplectic forms parametrized by (a, b) and the same Lagrangian on TS 2 containing two terms, of 
the form y 2 and (A(y x y),y x y), namely 

2T/(l+b/a) 2 =^ 2 Y 2 + (^(Yxi),Y><Y) • (3-22) 

For a/b = 1 (the case of equal spheres) we have /= 1 and (J.K) = 0, so we have the standard symplectic form in the base 
T*S 2 . Based on the classic integration by Jacobi of the ellipsoid geodesies problem, we could hope that the the Hamiltonian 
separates in sphero-conical coordinates X\,%2, so the Hamiton- Jacobi method would work. We were be lead to a disgusting 
failure. 
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3.5 Sphero-conical coordinates 

This section is also standard material for the experts. Sphero-conical coordinates appear naturally if one wants to express d3.22t 
in the form 2T = c\X 2 + C2^2 ■ To achieve this aim we need to diagonalize the positive quadratic form (Ax,x) , x 6 R 3 , A = 
diag(l\,l2,h) when restricted to subspaces defined by (x, y) = 0. This is a well known problem, going back to Lord Rayleigh. 
It boils down to extremize (Ax,x) subject to = 1 and (x,y) = 0. Let us apply the Lagrange multiplier method, with an 
ammended objective function 

U = \ L l A - HL*J - J ) -vLxjJj ■ (3-23) 

The first order condition VU X = yields 

Xj=fi - r * T , l<j<3. (3.24) 



Inserting in the constraints = 1 , L x ili = we get 



LA=0 (3-25) 



^VE^]. (3.26) 

This holds in n dimensions, but let us keep working with n = 3 for concreteness. If we order the parameters as I\ < I 2 < It,, 
equation B.25t has two roots, 

/, < Xi =X 1 (y;/i,/ 2 ,/ 3 ) <h 

I 2 < X 2 = X 2 (yJuh,h) <h 
(they can be found explicitly as the roots of a second degree equation in X thinking of the jj as parameters). 

For each choice of X, ( 13.25b represents an elliptical cone in the three dimensional space (Yi > Y2, 73)- It i s eas Y t0 see that the 
curves of intersection of these cones with the sphere (y,y) = 1 form an orthogonal system of coordinates. This is because we 
know that the two vectors J3.24t that solve d3.23t are perpendicular, and moreover, these vectors are manifestly parallel to the 
gradients of d3.25l l. Thus we can anticipate that 

(A(yxj),yxj)=b l X 2 1 +b 2 Xl . (3.27) 

In order to find these coefficients b\,b 2 , we need to express (Yi , y 2 , Ys) in terms of X\,X 2 . This is simply a matter of solving the 
linear system 

A(7?,Yi,YL} f = (l>0,0y 

with 

/ 1 1 1 \ 



A = 



(i,-x,)- x (h-Xj)- 1 (I 3 -Xj)- 1 



\ (i,-x 2 r x {i 2 -x 2 y l {i 3 -x 2 y l J 

Inverting A gives the desired explicit expression, 

{h-X^ih-Xa) {l2-X l )(I 2 -X 2 ) (I 3 -X l )(I 3 -X 2 ) 



\tu 72,73; 1 (I _ h] > (h _ h](h _ h] . (h _ h](I ._ h] j ■ v* 1 *) 



It follows that 



and 



{h-h){h-h) (h-h){h-h) {h~h){h-h) 
dy 1 Yi y 2 Y3 



1,2 (3.29) 



1 dXi 2 \{I2-Xi){l3-Xi) ' {h-Xi)(I 3 -Xi)' (h-XiXh-Xi) J ( - ' 

The remaining calculations can be done by brute force, but it was less painful to use a computer algebra system. The results are: 



Lemma 5 

(Ay,Y) = l/4-^ — ^ + 1/4- '7";. — XI (3.31) 



(X 2 -X])X] 12,1/4 {h-h)h j 2 

(Xi-Ii)(Xi-h)(U-I 3 ) 1 ' (h-Ij)(h-I 2 )(X 2 ~l3) 2 



m2 = 1/4 (^-/,)(t-/ 2 )(^-// ' + 1/4 ^-mC-m^-is) %l 
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Equation 0.3U indeed tell us that the ellipsoid geodesies problem separates (the standard metric x 2 in the ellipsoid £ jc? /Ij = 1 
is equivalent to the deformed metric (Ay,y) in the sphere Y.jj = 1 , Jj = xj/y/Tj). The standard metric J3.32b on the sphere, 
which gives a superintegrable system, also appears separable in conical coordinates (but in a somewhat cumbersome format). 

Let us now write the symplectic form d3 -21b using ^1,^2, Xi,Xq as coordinates. This requires expressing the following 
objects in terms of those: 



Lemma 6 



dS 2 = (y, iL x ^-)dXi AdX 2 = h Xl dXjdXg 

^1 a v v(h-h)(h-h)(h-h)(h-h)(h-hm-h) 



MM) = ^mrs fa {h %£l kl) + ^ + ^ {h -Sm-h l ) +cycUc: 



1 l(\i-ii)(\i-h)(h-h)^ i + l {h-h)(h-h)(h-h) ^ 



2V (Xi-/i)(/ 2 -A,l)(/3-Xl) 2 V (A.2-/i)(X2-/2)(/3-A. 2 ) 

/2 ^3 A hh 

yj+ti y?+Y3 y?+yj _ ^1+^2 

hh hh hh hhh 



F(y) = (l + (A* +^ 2 )(^2 +/* 2 ) (3.34) 



l-^+J lAa A l)| v( X l -/ l )(/ 2 -^ l )(/3-^l) l + (^2-/l)(^2-/2)(/ 3 -X2) 2 

Pl =3T/ ^ = ^ 1 + ^- X ') (X l -/ l )C?2-t)(/3-X l ) Xl ^ 

' 2 = 3r ^ 2 - ^ + ^-^ ^^-/o^-l)^-^) ^ (337) 

Putting everything together, we get: 

Proposition 10 Darboux coordinates in new time: 

®conformal = [(^1 +/* 2 ) (^2 + A* 2 )] ~ ®nh = (3-38) 

= [(Xi+pb 2 )(X 2 +pb 2 )}^ (dpidkt + rf P2 ^2 + \ ( 1 - fe/a) ( ^ - ^ )dh dX 

= d(P x dXi +P 2 dX 2 ) Darboux! (3.39) 
Proposition 11 The Hamiltonian in spheroconical coordinates: in new time 

dx/dt = , F(j) = (1 + -J 2 (Xi +pb 2 )(X 2 +/Jb 2 ) (3.40) 

defining 

Pi = [(Xi+pb 2 )(X 2 +pb 2 )}^ Pl , P 2 =[(X l +pb 2 )(X 2 +pb 2 )}^ P2 . (3.41) 

the Hamiltonian becomes 

„2 2 p2 p2 

2H = — + — = ! ~ — + 1 rr— (3-42) 

c i c 2 [(X l +pb 2 ){X 2 +^b 2 )]-c l [(X 1+ pb 2 )(X 2 +vb 2 )]-c 2 

with 

Cl = 7(1+-) (A-2 - Aj)-^ — — r— —,c 2 = - 1 + -) (A2-A.1 -r — — — — jpr. (3.43) 

4 a (/.I (I 2 — X\ ) (/} — Xi) 4 a (X 2 - h)(X 2 - I 2 )(h -X 2 ) 
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4 Some directions for research 
4.1 Nonintegrability. 

If the rolling ball has its mass concentrated at the geometric center, then the inertias /; vanish. Hence, the problem with Ij = zlj, 
is a perturbation of the geodesies on the sphere. As a possible attempt to prove nonintegrability rigorously is to compute 
the monodromy of the linearized equations along a great circle (solution of the unperturbed problem) and apply Melnikov or 
Moralis-Ramis methods. Parhaps this may be done directly in y 6 R 3 . In sphero-conical coordinates, using rescaled variables 
Xj = sXj and new time with factor 

F(y) = (1 + -f (eXi +iib 2 ){zk 2 + M b 2 ) (4.1) 
a 

we obtain (dropping the bars) the following Hamiltonian perturbation of the geodesies on the sphere: 

p2 p2 

2H = ! = 1 ? = (4.2) 

[(eX,i + A* 2 ) (e^ 2 + vb 2 )] - fli [(Ai + a* 2 ) (eX 2 + ^ 2 )] ^ a 2 



with 



ai = J 1 + 7l> ^- Xl) (X I -/ 1 )(/ 2 -3L I )(^-3L 1 )' fl2= 4( 1 + ^ (X2 - Xl) (X 2 -/ 1 )(^-/2)(/3-X 2) - (43) 



4.2 Rolling with sliping but no twisting and Rattleback phenomena 

From the mathematical point of view, it is legitimate to study systems subject to the no-twist constraint, sliding allowed. 
The question is whether these constraints could be (even approximately) realized in practice. Nonholonomic constraints are 
approximately realized as limits of certain collisional systems. It would be interesting to conceive a "peg-leg" version for 
no-sliping and the no-twist constraint^. 

We believe that the dynamical torque l |2.25t the heart of the rattleback rotation reversal phenomenon. There is a extensive 
literature on the topic, but we would like to see an explanation based on the joint geometries of Z 2 , say with E] being the plane. 
Is there a rubber rattleback? Here the idea is not looking at rotations (which are forbidden) but at translation reversal. 



4.3 Skiding dynamics: mathematical hockey 

As far as we know holonomic systems T(q\,R,q\,R) given by l !2.4| >, describing the motion of a body <B skiding over a surface 
Zi have not been studied systematically. In general the rotational and surface dynamics should be coupled and some interesting 
effects of mathematical hockey may be anticipated^!. A first project in this direction is to explore symmetries (as we discussed 
briefly). 



4.4 Rubber rolling of surfaces of revolution over a sphere 

Conditions under which a generalized Chaplygin system has an invariant measure are given in 1 171 . It would be interesting to 
make explicit which condition results for marble rolling of an arbitrary surface E 2 over a sphere or a plane. Also, implementating 
proposition [3] on exemples such as bodies of revolution and ellipsoids is in order, as well as the higher dimensional extensions 
13011361 . These situations have been studied in the marble rolling context, see tables 1 and 2 of 1 10 1. On another talk, rubber 
rolling of an homogenous sphere over arbitrary surfaces Zi , see 1111 for the marble counterpart. 



4.5 Suslov rolling 

While the usual no-slip constraint is enforced, the no-twist constraint CO3 = is replaced by il^ = 0, that is, one imposes that the 
third component of angular velocity in the body frame should vanish. This constraint can be implemented by a inner spherical 
support. 



4.6 Rubber rolling and Yang-Mills 

Among all convex bodies <S with the same area, rolling over a sphere, which one minimizes the Yang-Mills functional? This 
functional, important for gauge theories, is given by 




'We thank Andy Ruina for a discussion on this subject; see http://ruina.tam.cornell.edu/ 

Hockey (see http://www.exploratorium.edu/hockey/ 1 could be even more exciting with a strictly convex puck. 



REFERENCES 



20 



where K is the curvature of the connection. Here we restric ourselves to rubber connections ty 3 on the bundle 50(3) x S 2 . A 
nice exercise is to compute the rubber connection 1-form and its curvature, for a tri-axial ellipsoid rolling over a sphere. 

A reasonable conjecture is that the sphere minimizes the YM functional. Clearly, the minimum attains the value zero when 
the area of B is the same as the area of S a . 

Yang Mills connections over Riemann surfaces were considered on an important paper by Atiyah and Bott \ 3\. Particularize 
to rolling problems when the base is a surface of constant negative curvature. 
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